CMB Optical Depth Measurements: Past, 

Present, Future 



Brian Keating^ and Nathan Miller^ 

:-\ ^Department of Physics, University of California, San Diego, Mail Code O424, La 
^ ■ JoUa, CA, 92093-0424, United States of America 

. E-mail: hkeating@ucsd.edu, nmiller@physics.ucsd.edu 

(N 
o 

CD 

Q 

00 



> 



Abstract 



The polarization of the cosmic microwave background (CMB) is encoded with ex- 
actly the same cosmic information as the CMB's temperature anisotropy. However, 
Q"^ ' polarization has the additional promise of accurately probing the reionization his- 

. tory of the universe and potentially constraining, or detecting, the primordial back- 

i ground of gravitational waves produced by inflation. We demonstrate that these two 

' CMB polarization goals are mutually compatible. A polarimeter optimized to detect 

■ the inflationary gravitational wave background signature in the polarization of the 

CMB is well situated to detect the signatures of realistic first-light scenarios. We 



also discuss current results and prospects for future CMB polarization experiments. 
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^ ' 1 Introduction 

The CMB possesses a wealth of cosmological information owing to its origin 
400,000 years after the Big Bang - an epoch when linear cosmological per- 
turbation theory holds. However, the CMB also is imprinted by "first-light" 
sources at later times that reveal the dynamics of reionization. Even detect- 
ing the transition from partial to full-reionization in (physically more plau- 
sible) non-instantaneous reionization scenarios is conceivable. The first-year 
WMAP results demonstrate a nearly 5a detection of a reionization feature 
in the temperature-polarization cross-correlation power spectrum. This detec- 
tion, along with WMAP's detection of the Thomson optical depth, was the 
first demonstration of CMB polarization's ability to go "beyond parameter 
estimation" - making precision measurements of a cosmological epoch which 
was previously shrouded in mystery. WMAP's detection was also a tantalizing 
hint at CMB polarization's promise to probe vitally important cosmological 
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parameters that are essentially unobservable using temperature anisotropy 
alone. 

However, as with many groundbreaking experiments, WMAP's results have 
generated even more intriguing questions. In particular, why is the Thomson 
optical depth so much higher for the CMB than that implied by the QSO 
Gunn-Peterson measurements, e.g. [1,2,3]? Ultimately, the optical depth due 
to first light sources measured by CMB polarimeters must be reconciled with 
the Gunn-Peterson absorption feature measurements in the spectra of distant 
quasars. Several attempts at this reconciliation have been proposed, some of 
which feature multiple reionizations [4], or a period of partial ionization [5]. 
Interestingly, it appears that these reionization scenarios, which translate to 
structure formation scenarios, can indeed be constrained observationally by 
GMB polarization. 

In this paper we discuss the prospects of upcoming ground, balloon, and space 
CMB polarimeters to measure both the Thomson optical depth and constrain 
new features in the polarization power spectra characteristic of realistic reion- 
ization scenarios. The implications of intrinsic systematic effects such as fore- 
grounds and cosmic variance are also discussed. 



2 CMB Polarization 



Much of the wealth of cosmological information encoded in CMB polarization 
is only observable at medium to large angular scales (> 0.5°). Essentially all 
information on the classical cosmological parameters of the ACDM model, 
e.g., the baryon density, Qh, the Hubble parameter, h, the scalar power spec- 
trum's power law index, Hs, the amplitude of the primordial power spectrum, 
A, the tensor-to-scalar ratio, r, etc, can be obtained from CMB polariza- 
tion at angular scales > 0.5°. Because the temperature anisotropy is much 
larger than the polarization signal it can produce higher confidence measure- 
ments of the classical cosmological parameters than polarization. However, 
there arc several parameters which only polarization can reveal. In particular, 
only CMB polarization can probe the imprint of the infiationary gravitational 
wave background (GWB) on the CMB. This polarization signature occurs at 
angular scales of ~ 2°. Reionization primarily affects CMB polarization on 
angular scales > 4°. We note that the grad-mode polarization, Cf, peaks at 
i ~ 1000 corresponding to ~ 10' scales, and that gravitational lensing of the 
CMB grad-mode polarization by large scale structure produces non-negligible 
curl-mode polarization Cf at the same ~ 10' scales [6] . 
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3 Current Results 
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Fig. 1. Measurements of the large-angular scale gradient-mode polarization power 
spectrum Cf for 0.5° <^ £ < 400. The solid line is the polarization power spectrum 
for the WMAP best fit cosmological model, with r = 0.17 [7]. 

Figures 1 and 2 show the current detections of the grad-mode polarization and 
the polarization-temperature cross correlation. Currently, only DASI[8] and 
BOOMERANG [9] have detected grad-mode polarization at angular scales < 
0.5° i < 400, and no experiments have detected the grad-mode polarization 
for i < 100. The only other experiments that have detected the grad-mode 
polarization are CBI[10] and CAPMAP[11], though at smaller angular scales 
(higher i) than are considered in this paper. 

The primary effects of reionization are encoded in the grad-mode and tem- 
perature polarization cross-correlation at £ < 50. In the absence of detections, 
the most stringent constraint in this range of multipoles is currently Cf < 8 
fiK at 95% confidence reported by P0LAR[12] for 2 < £ < 20 (assuming no 
B-modes). There are many more detections of the temperature-polarization 
cross correlation than of the grad-mode polarization as figure 2 demonstrates. 
WMAP [13] reports a large number of highly significant detections, especially 
at low-i due to WMAP's ability to map most of the sky. CBI, DASI, and 
BOOMERANG [14] have also detected the cross-correlation spectrum, mainly 
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Fig. 2. Measurements of the temperature-polarization cross-correlation power spec- 
trum Cj^ . The solid line is the power spectrum for the WMAP best fit cosmological 
model, with r = 0.17 [7]. 

at much smaller angular scales than WMAP. A complete description of reion- 
ization will require significant detections of CMB E-mode polarization for 
£ < 50. An ancillary benefit of reionization is that it boosts the primary curl- 
mode power spectrum significantly near I = 10. Due to reionization, a more 
stringent limit on the tensor-to-scalar ratio, r, in the presence of lensing can 
be obtained than that calculated in [15,16]. 



4 The effect of reionization 



Reionization produces free-electrons which Thomson-scatter CMB photons 
and produce CMB polarization. The scattering damps the CMB temper- 
ature in direction fi as T'(n) = e~^T(n), where r is the Thomson opti- 
cal depth. Reionization, therefore, damps the anisotropy power spectrum as 
Cj = e~^'^Cj. The effect of r on the temperature anisotropy power spectrum 
is completely degenerate with the primordial power spectrum's amplitude, 
A. Fortunately, a new feature, e.g. [17,18], in the polarization power spectra 
occurs at large angular scales which breaks the degeneracy between normaliza- 
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Fig. 3. Temperature (color scale) and polarization (lines) realizations in real-space 
smoothed with a 4° beam. The left panel shows a simulation with no reionization 
and the right panel shows the effects of instantaneous reionization with r = 0.17. 
Reionization does not noticeably affect the large scale temperature pattern but it 
produces roughly ten-times larger polarization at large scales. These figures have 
no curl-mode power. The maximum gradient-mode polarization in the right panel 
is ~ 800 nK. Figures courtesy of Eric Hivon. 

tion and optical depth. The degeneracy (for the CMB temperature anisotropy) 
and its breaking (for CMB polarization) is nicely illustrated in (real-space) in 
figure 3. 



5 A More Realistic Cosmic Variance Limit. 

An all-sky experiment must confront galactic contamination, which can only 
be subtracted to finite precision. For example, the WMAP team used multiple 
frequency channels in combination with a sky-cut that left 85% of the sky re- 
maining [13]. This cut also introduced correlations between multipoles at the 
few percent-level, and a total anti-correlation of 12.4%. Even with sophisti- 
cated foreground modelling and multiple frequency coverage the challenges of 
achieving polarization fidelity at the 100 nK level over > 10° angular scales - 
required, for example, to detect the signatures of non-instantaneous reioniza- 
tion scenarios (e.g, [5,19]) - are significant. 

As a simple, but admittedly conservative, foreground mitigation technique we 
consider a full-sky map and and omit galactic latitudes within |6| < bcut of the 
galactic plane. This leaves two circular "cap" regions remaining. This strategy 
affects the recovery of reionization information in two ways. First it reduces 
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the resolution in £-space compared to a true full-sky experiment, leaving a 
residual correlation between multipoles. Secondly, it increases cosmic variance 
on the measurement of Cf , since the amount of observed sky, fgky decreases: 
5Ci oc 

J sky 

The effect of this simple foreground mitigation strategy is apparent in parame- 
ter estimation. If one (conservatively) insists on using uncorrelated multipoles 

^cutoff 

to estimate parameters, the integrated signal-to-noise ratio ^ TT^e 

£=2,4,6... 

duced significantly. Here 5Cf is the cosmic variance limited uncertainty on 
Cf, including the inflation of the cosmic variance due to the cut itself, and 
^cutoff is the highest multipolc afi'orded by the angular resolution of the tele- 
scope or other filtering. 

As a numerical example, we consider a galaxy-cut which excludes data with 
bcut = 20°, corresponding to fsky = 75% - not far off from WMAP's fsky 
However, since we insist that the total anti-correlation between multipoles is 
<^ 1%, the ^-space resolution degrades from Si — 1 to S£ — 2. The quadrupole, 
and all even multipoles are completely uncorrelated. 

We find that the signal-to-noise is more than halved relative to a true all-sky 

^cutoff 

experiment ^ — ^ . The halved signal-to-noise ratio of the power spectra 

£=2,3,4... 

does not necessarily translate into a doubling of parameter errors, however. A 
Fisher matrix analysis, e.g. [20], must be performed to assess the degradation 
of parameter sensitivity. 

A parameter estimation indicates that the new cosmic variance limit on r 
from the hcut — 20° cut-sky data is a-,- — 0.0048. This is compared to our 
Fisher matrix results for the ideal full-sky cosmic variance limit: a-r — 0.0029 
- an increase of 65%. This latter limit agrees nearly exactly with the results 
of Holder et al. (2003) [19]. Due to the nearly degenerate correlation between 
r and the amplitude of the primordial power spectrum, the increase in 
increases the uncertainty on A as well. Figure 4 shows the individual and 
joint hkelihoods for r and A for the full-sky and bcut — 20° simulations. 

Additionally, more realistic cut-sky observing strategies have lower £-space 
resolution than all-sky surveys which may preclude the detection of features in 
the polarization power spectra that characterize the dynamics of reionization 
[5,19] . This may ultimately limit the precision with which the dynamics of non- 
instantaneous reionization can be reconstructed. Furthermore Kaplinghat et 
al.[5] estimate that the sensitivity of a CMB polarimeter to the redshift of 
the onset of reionization is cr^^^ = 333aT-/y^. This implies that cr^^. for a 
more realistic observation strategy (i.e., cut-sky) will be determined 65% less 
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Fig. 4. Individual and joint likelihood functions for the cut-sky and full-sky obser- 
vations discussed in the text. The errors on r and A from the ID likelihoods for 
the bcut = 20° simulation (solid) are inflated by ~ 65% over the true full-sky ex- 
periment (dashed) due to increased cosmic variance and decreased multipole space 
resolution. The 2D joint likelihood contours show the full-sky survey's error volume 
(white) inside the the cut-sky data's region (black). The contours indicate where 
the likelihood falls to l/s/e of the maximum likelihood for each experiment. 

precisely than an ideal truly all-sky measurement. These effects will thus be 
more pronounced as the sensitivity of CMB polarimeters increase. 

This more realistic cosmic variance limit may also strengthen the motivation 
for ground-based CMB polarization probes of reionization which, due to their 
significantly reduced sky coverage relative to a satellite or balloon, have re- 
ceived comparatively little attention. 



6 A Ground Based Reionization Probe 



In this section we explore how well a ground based polarimeter can probe 
reionization. We seek to quantify the observational requirements to detect r, as 
well as to constrain realistic reionization models consistent with both WMAP 
and quasar absorption measurements. In particular we assess the requirements 
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to detect a two-step reionization model suggested by [5] . 

Generally speaking, an instrument optimized to detect gravitational waves will 
significantly constrain the detailed reionization history, whereas the converse 
is not necessarily true. For fixed system sensitivity this has implications for 
both the optical design and observation strategy of the experiment. 

Here we will consider a medium angular-scale, ground-based polarimeter's 
ability to detect r and constrain non-instantaneous reionization scenarios. As 
a toy experiment, we consider a polarimeter with 1° resolution and system 
sensitivity (NETgygt = '^^"''pixel/\/^^^pixd ^ 70^Ks^/'^, observing for one year. 
The system sensitivity for polarization is ^/2 times higher NEQ = ^/2NET. 
These detector requirements are well within reach of the current generation 
of CMB polarimeters. 

The problem of optimizing of CMB polarimeters for detection of the infla- 
tionary GWB was first studied by [21], and polarimeter optimization in the 
presence of the gravitational lensing foreground has been studied by Lewis, 
Challinor, and Turok [22] . Lewis et al. find there is a broad maximum sensitiv- 
ity achieved by mapping a circular cap of radius © = 20°, which corresponds 
to an ^-space resolution of 5£ ~ 5. In the absence of galactic foregrounds, 
our toy experiment is capable of detecting the curl mode signal imprinted by 
primordial gravitational waves with a tensor-to-scalar ratio r = 0.12 at 95% 
confidence, with no priors. This is roughly ten-times lower than the limit on 
r from WMAP, also with no priors [7]. This detection (or hmit) would per- 
mit discrimination between a broad class of inflationary models [23,24,25]. 
An exciting ancillary beneflt of this experiment is that it is also sensitive to 
reionization. 

To assess the power of the toy experiment to constrain the Thomson optical 
depth we consider the susceptibility of the E-mode power spectrum to changes 
in r: dCf /dr. From the grad-mode (Fig. 5) and curl- mode derivatives (Fig. 
6) we see that essentially all the information on r comes from i < 50. If we 
note that the peak in the curl-mode power spectrum occurs at £ ~ 90, we 
see that the scale for completely resolving the effects of r is only two-times 
larger than that required to resolve the curl-mode structure. While most of 
the r-constraint comes from i < 10, a non-negligible amount comes from 
10 < £ < 40. There is significant power beyond the peak of the low-£ bump. 
A Fisher matrix analysis for the toy model indicates its sensitivity to r is 
(7,- = 0.08, competitive with the one-year WMAP results [7] with no priors. 

It is instructive to explore how the sensitivity to r chang function of 

sky coverage. In table 1 and figure 7 we demonstrate how r— sensitivity, o",-, 
depends on £-space resolution. The £-space resolution scales as Si ~ 18O°/20. 
With one exception, for all of the observation strategies we consider we ex- 
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Fig. 5. The susceptibility of the E-mode power spectrum to changes in r. The 
derivatives were computed using a second-order finite difference with steps in r 
of 5t = 0.001, around a fiducial value of r = 0.17, consistent with WMAP. The 
dominant effect of r is at ^ < 10, but a non-trivial amount remains for 10 < £ < 40 
where cosmic variance effects can be dramatically reduced, even for a ground based 
experiment. 




Fig. 6. Same as Fig. 5, but for curl-mode polarization. 



elude the first i-hin which is, in principle, observable. This results in V^i- 
hin— 251. This conservatism is frequently justified since the first ^-bin often 
is contaminated by an imperfect foreground template and/or by subtraction 
of the DC-level of the map. WMAP's purportedly "anomalously-low" first l- 
bin, ^ = 2 has received much attention, perhaps prematurely. However, some 
have speculated on the potential for foreground contamination [26] of, or even 
intriguing new physics[27] in, the quadrupole's low value. 
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Table 1 

Sensitivity to optical depth, r, versus observation strategy, parameterized by radius 
of circular cap, 0, observed by the toy experiment described with system NET^j^st = 
TO/iKs-*^/^. For all of the observation strategies, except 9 = 9°, the first ^-bin is 251. 




Fig. 7. Signal-to-noise ratio for a ground based polarimeter to detect rcionization, 
assuming r = 0.17. The plot shows the effect of increasing the angular radius 
of the circular survey, 6, with constant system sensitivity NET^jy^j = 70/uKs^/^. 
All strategies except © = 9° exclude the first multipole bin from the parameter 
estimation. Lewis, Challinor, and Turok[22] show that G = 18° is close to optimal 
sky coverage to detect the inflationary GWB using CMB curl-mode polarization. 

COSMOMC is a Monte Carlo approach [28] useful for evaluating sensitivity to 
cosmological parameters. Figure 9 shows the parameter estimation results of 
COSMOMC for the toy experiment. Due to the cut-sky and the intrinsically 
low signal-to-noise ratio of reahstic reionization probes, such as WMAP, it is 
unlikely that the marginalized likelihood functions for r will be convincingly 
gaussian in nature. Indeed, Spergel et al. explicitly points out [7] that r is the 
only parameter with a highly non-gaussian form. COSMOMC facilitates the 
exploration of the details of the likelihood curves and the effects of adding 
priors on the Thomson optical depth. 
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Fig. 8. The effects of including prior information on r on the hkehhood functions for 
the toy experiment described in the text. The Ukehhoods were computed with sim- 
ulated data for the toy experiment, using priors ranging from r > 0.02 to r > 0.07. 
The likelihood curve for WMAP with no priors is shown for comparison as the 
dashed curve. The upper la range of the likelihood curve is robust to changing the 
prior, whereas the lower range is more sensitive. 

Figure 8 shows the hkelihood functions for r for the toy experiment with 
G = 18° as a function of r-priors ranging from < r < 1 to 0.07 < r < 1 
(weaker constraint than the SDSS Gunn- Peterson observations provide). The 
toy-experiment is also compared to the WMAP-only likelihood result in the 
same figure. 

As discussed above, Kaplinghat et al. demonstrate that the ability of a CMB 
polarimeter to constrain the redshift of reionization is cr^^. 
Thus, conservatively, if reionization is consistent with the WMAP la lower 
limit [zri — 13), the toy experiment with (0 = 18°) could establish (with 
> 95% confidence) that reionization was earlier than z = 6.28. This is the 
redshift of the SDSS quasar neutral hydrogen fraction determination from 
Gunn-Peterson QSO absorption features. We note that these results are for 
analysis of the grad-mode data only. The toy experiment can either use its own 
data to measure the GMB temperature anisotropy or WMAP's. Including the 
temperature-polarization cross- correlation data the precision on r improves 
by ~ 20%. Finally, since r and the primordial power spectrum amplitude A 
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Fig. 9. Individual and joint likelihood functions for the toy experiment discussed 
in the text, with fsky = 2.4%. The uncertainties from the marginalized likelihoods 
are r = 0.12 ± 0.038 and A = 0.81 ± 0.15, given a prior that r > 0.04. The joint 
confidence intervals correspond to l/-v/e (inner) and 1/e (outer) of the maximum 
likelihood. This experiment can determine the redshift of the onset of reionization 
with a precision oi 5z < 3 if reionization occurred at z > 10. 

are exactly degenerate, = 2(7,-, and the toy experiment can achieve 6% 
fractional uncertainty on A. 



7 Prospects for the future 

WMAP's on-going measurements of CMB polarization are expected to con- 
tinue until at least 2006. In 2007 Planck will be launched and will start making 
full-sky observations of the polarization of the CMB. Planck is expected to 
achieve near-cosmic variance limit on r. NASA's Beyond Einstein initiative 
features a CMB polarimeter called the Inflation Probe [29]. This instrument 
is designed to detect the signature of inflationary gravitational waves over 
a wide range of inflationary energy scales. This experiment will achieve cos- 
mic variance limited precision on r, which is essentially the same as Planck's 
sensitivity [30]. The European Space Agency plans a large angular scale po- 
larimeter deployed on the International Space Station called SPORT, which 
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would also be sensitive to reionization[31]. 



8 Discussion 

This paper has emphasized that there is more to CMB polarization than sim- 
ply the lifting of cosmic degeneracies. Upcoming experiments will be able to 
probe the details of reionization and the signature of the gravitational wave 
background simultaneously. We have emphasized that there is far more in- 
formation in CMB polarization data than whether or not two power spectra 
can be distinguished from one another. However, conservative foreground mit- 
igation techniques appear to degrade the precision with which the epoch of 
reionization, and even the Thomson optical depth itself, can be obtained by a 
cosmic variance limited experiment. Whether such conservatism is warranted 
is debateable. But when probing signals below lOOnK level it is clear that 
foreground mitigation is vital. More sophisticated techniques may restore lim- 
its closer to the full-sky cosmic variance limit. Detailed simulations are war- 
ranted. A current generation ground based polarimeter is capable of probing 
the details of the reionization epoch as well as the signature of inflationary 
gravitational waves. Clearly, high-fidelity CMB polarization observations are 
poised to enhance and extend our knowledge of the cosmic Dark Ages. 
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